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Abstract 


The effect of gluon shadowing on charm quark production in large im¬ 
pact parameter ultrarelativistic heavy ion collisions is investigated. Charm 
production cross sections are calculated for a range of non-central impact pa¬ 
rameters which can be accurately inferred from the global transverse energy 
distribution. We show that charm production is a good probe of the local 
parton density which determines the effectiveness of shadowing. The spatial 
dependence of shadowing can only be studied in heavy ion collisions. 
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I. INTRODUCTION 


Deep inelastic scattering experiments using nuclear targets showed that the quark and 
antiquark distribution functions are modified in the nuclear environment [1] and hence are 
different in heavy nuclei than in free protons. It is not unreasonable to expect the nuclear 
gluon distributions to be affected at least as much as the quark distributions. However, little 
is known about the nuclear gluon distribution because the gluon distributions can only be 
indirectly probed. Gluon-dominated production processes, such as J/ij) and heavy quark 
production, can provide an indirect measure of the nuclear gluon distribution. Since the 
J/i[} is more strongly affected by absorption processes than charm quarks, evident from their 
respective A dependencies [|],[|, charmed quark production provides a cleaner determination 
of the nuclear gluon distribution. 

To date, all measurements and indirect determinations of nuclear parton distributions 
have been insensitive to the position of the interacting parton within the nucleus. However, 
there is no reason to expect the parton momentum distributions to be constant within the 
nucleus. They should at least vary with the local nuclear density. If shadowing is due to 
gluon recombination, the position dependence could be quite strong f|. One way to probe 
the position dependence of the shadowing is to measure cc production over a wide range 
of impact parameters, thus scanning gluon localization in the nucleus. The charm rate has 
been shown to be large in central collisions ||, here we will show that these studies are also 
feasible at large impact parameters. 

This paper thus proposes a method for measuring the position dependence of the gluon 
momentum distribution in heavy nuclei. We show that the charm production rates in non¬ 
central collisions are sensitive to the details of the gluon distribution and its position de¬ 
pendence. We use two different parameterizations of nuclear shadowing along with two 
parameterizations of the position dependence of the shadowing to calculate charm produc¬ 
tion in 100 GeV per nucleon Au+Au collisions at the Relativistic Heavy Ion Collider (RHIC) 
[6|, now under construction at Brookhaven National Laboratory. However, the techniques 
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discussed here should also be applicable to cc and bb production in Pb+Pb collisions at 
the CERN Large Hadron Collider (LHC). The charm quark production rate and px spectra 
are calculated as a function of impact parameter, b , for non-central collisions with impact 
parameters greater than the nuclear radius, Ra- 

For this study, we need to select events according to impact parameter. Although the 
impact parameter of the collision is not directly measurable, it may be inferred from the total 
transverse energy, E T , of the event 0. We discuss the relationship between E T and b and 
present calculations showing that, for a given Et, the impact parameter can be measured 
relatively accurately. Additional input, such as a measurement of nuclear breakup, through 
the use of a zero degree calorimeter, can refine this estimate. 

Section 2 summarizes the calculations of cc production in peripheral collisions including 
a discussion of the nuclear parton shadowing and its possible spatial dependence. Section 
3 discusses the relationship between transverse energy and impact parameter. Section 4 
presents the numerical results for the charm production rates and px spectra for two ranges 
of non-central impact parameters. We demonstrate how these rates are sensitive to the 
nuclear gluon distribution. Our results are put into an experimental perspective in Section 
5. Finally, Section 6 draws some conclusions. 


II. CC PRODUCTION 


To study the effects of shadowing on cc production in peripheral collisions, we emphasize 
the modifications of the parton distribution functions due to shadowing as well as the location 
of the interacting parton in the nucleus. We discuss the method used to calculate cc pair 
production and introduce two parameterizations of nuclear shadowing. We also describe two 
models of the spatial dependence of the shadowing. 

The double differential cross section for cc pair production by nuclei A and B is 


E c Ec 


d(JAB 

d 3 p c d 3 Pcd 2 bd 2 r 


( 1 ) 
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Y2 J dz dz' dx i dx 2 F 1 A (x 1 , Q 2 , r, z)Ff(x 2 , Q 2 , 6 

i,3 J 


r, z')E c E-, 


da ij (x 1 P 1 ,x 2 P 2 ,rn c , Q 2 
d 3 p c d 3 pc 


Here i and j are the interacting partons in the nucleus and the functions F t are the number 
densities of gluons, light quarks and antiquarks evaluated at momentum fraction x, scale 
Q 2 , and location r, z. (Note that f is two-dimensional.) The short-distance cross section, 
a.ij, is calculable as a perturbation series in a s (Q 2 ). 

At leading order (LO), cc production proceeds by two basic processes, 


q + q —> c + c (2) 

g + g ^ c + c . (3) 


The LO cross section, 0(a 2 ), can be written as 
dcTAB 


d 3 p c d 3 pcd 2 bd 2 r 


( 4 ) 


— dz dz' dx i dx 2 C(xi, x 2 , Q 2 , r, z, b — r, z') 5 4 : (xiPi + x 2 P 2 — p c — pz) 
J 2 tt 


where y/s, the parton-parton center-of-mass energy, is related to y/S, the hadron-hadron 
center-of-mass energy, by s = XiX 2 S > 4 m 2 , where the momentum fractions, X\ and x 2 , are 


xi ,2 = ~^(e ±y + e ±y ) , (5) 

V s 

and ttit = \]rn 2 + ~jFr- Tl ie target fraction, x 2 , decreases with rapidity while the projectile 
fraction, xi, increases. Here, the intrinsic transverse momenta of the incoming partons has 
been neglected. The convolution of the subprocess cross sections with the parton number 
densities is contained in C(x i, x 2 , Q 2 , r, z,b — r, z') where 


C(xi,x 2 , Q 2 , r, z, b - r, z') = 


( 6 ) 


Yll F f( x i> Q 2 ’ F z)F?(x 2 , Q 2 , b - r, z') + F- A (x i, Q 2 , r, z)Ff(x 2 , Q 2 , b - r, z')\ 
<? 

Q 2 , F z)Ff{x 2 , Q 2 ,b~F ■ 


r\ I d(Jqq 


dt 


Four-momentum conservation leads to the rather simple expression 
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d&AB 


= f dz dz' x 1 x 2 C{x 1 ,x 2l Q 2 ,f, z,b - f, z') , 


(7) 


dp^dydyd?bd 2 r 

The LO subprocess cross sections for cc production by qq annihilation and gg fusion, ex¬ 
pressed as a function of uit, y, and y, are 0 


ddgg 7ra 2 cosh (y — y) + m 2 /m. 


dt 


9 raj, (1 + cosh (y — y)) 3 


dd g g 7T a 2 8 cosh(y — y) 


1 ( u _x , 2 ml 
coshfy — y) + 


2 mi 


rrin 


( 8 ) 


(9) 


dt 96 r?4 (1 + cosh(y — y)) 3 y " ' m ^ ,,o T 

Leading order calculations tend to underestimate the measured charm production cross 
section by a constant factor, usually called a K factor, 

^exp (AB > Cc) 


/C L ° = 


( 10 ) 


exp a LO (AB ^ cc) ’ 

The next-to-leading order (NLO) corrections to the LO cross section have been calculated 
TOlJTHl and an analogous theoretical K factor /b th can be dehned from the ratio of the NLO 


to the LO cross sections, 


K t h = 


cr NLO (AB —> cc) 


( 11 ) 


cr LO (AB —>■ cc) 

where ctnlo is the sum of the LO cross section and the 0(a s ) corrections. 

Previously JT2|] , the NLO calculations were compared to the cc total production cross 
section data to fix m c and Q so that K^° ~ 1 to provide a more reliable estimate for 
nuclear collider energies. Reasonable agreement with the measured total cross section was 
found for m c = 1.2 GeV, Q = 2 m c for MRS D—' |L3| and m c = 1.3 GeV, Q = m c for GRV 
HO [|Id|1 . We choose different scales for the two sets[] because of the different initial scales of 
the two parton distributions. The MRS D—' distributions have Qq M rs = 5 GeV 2 ; we choose 
Q = 2 m c so that Q 2 > Qo.mrs- The GRV HO sea quark and gluon distributions are valence¬ 
like at low x and <3 o,grv = 0-3 GeV 2 . We can then use Q = m c because m 2 > Qo,grv- 
However, below Q 2 ~ 5 GeV 2 the gluon distribution is still somewhat valence-like. 


J These structure functions can be found in the CERN PDFLIB 15]. 
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When calculating inclusive distributions rather than total cross sections, it is more ap¬ 
propriate to choose Q oc m^, particularly when p T > m c since a constant scale introduces 
unregulated collinear divergences B- Therefore, we take Q = 2 rriT for the MRS D— 7 dis¬ 
tributions and Q = tut for the GRV ffO distributions. Both sets of parton densities result 
in a NLO total cc production cross section of ~ 350 pb in pp collisions at a/s = 200 GeV. 

The differential Ji t h for the charm quark p B distribution, the pair mass distribution, and 


the charm quark and cc pair rapidity distributions are nearly constant at RHIC energy |]T6 


They are also essentially independent of the parton density. The value of Jith is determined 
by a comparison of the NLO and LO total cross sections. Our LO calculations, eq. (4), are 
multiplied by the appropriate K th found for the specific parton density: 2.5 for the MRS 
D—' distributions and 2.9 for the GRV HO distributions. 

The nucleon parton densities are only a part of the space-dependent nuclear number 
densities, F A (x, Q 2 , r, z ), introduced in eq. (1). We have assumed that these nuclear number 
densities factorize into nuclear density distributions, independent of x and Q 2 , the nucleon 
parton densities, independent of spatial position and A, and a shadowing function that 
parameterizes the modifications of the nucleon parton densities in the nucleus, dependent 
on x, Q 2 , A and location, 

F a (x, Q\ f, z) = p A (s)S i (A , x, Q 2 , f, z)f[(x, Q 2 ) (12) 

Tf (x,Q 2 ,6- r,V) = p B (s')S l (B,x,Q 2 ,b-r,z')f?(x,Q 2 ) , 

where s = a /r 2 + z 2 , s' = \J\b — r| 2 + z 1 2 and ff are the nucleon parton densities. We 
assume that z and z 1 are uncorrelated. The collision geometry in the plane transverse to 
the beam is shown in Fig. 1. 

A three parameter Wood-Saxon shape is used to describe the nuclear density distribu¬ 
tion, 

1 + lu(s/R a ) 2 


Pa(s) = po¬ 


lls) 


1 + exp((s - R A )/d ) 

where R A is the nuclear radius, d is the surface thickness, and u allows for central irregu¬ 


larities. The electron scattering data of Ref. [ 17] is used for R A , d, and to assuming that the 
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charge and matter density distributions are identical. The central density, po, is found from 
the normalization / d 2, rdzpA{s ) = A. For gold, u = 0, d = 0.535 fin, R A = 6.38 fm, and 
po = 0.1693 fm- 3 . 

If the parton densities in the nucleon and in the nucleus are the same, then 
S l (A,x,Q 2 ,r, z) = 1. We will use this as a baseline against which to compare our results 
with shadowing included. 

We now discuss our choices of the shadowing parameterizations used in our calculations, 
independent of the position. Measurements of the nuclear charged parton distributions by 
deep-inelastic scattering on a nuclear target and a deuterium target, show that the ratio 
Rf 2 = F 2 4 / F[ j has a characteristic shape as a function of x. The region below x ~ 0.1 
is referred to as the shadowing region and the region 0.3 < x < 0.7 is known as the EMC 
region. In both regions a depletion is observed in the heavy nucleus relative to deuterium and 
Rf 2 < 1. At very low x, x ~ 0.001, Rf 2 appears to saturate^ [[HJ. Between the shadowing 
and EMC regions, an enhancement, antishadowing, occurs where Rp 2 > 1. There is also an 
enhancement as x —> 1, assumed to be due to Fermi motion of the nucleons. The general 
behavior of Rp 2 as a function of x is often referred to as shadowing. Although this behavior 
is not well understood for all x, the shadowing effect can be modeled by an A dependent 
fit to the nuclear deep-inelastic scattering data and implemented by a modiheation of the 
parton distributions in the proton. We use two different models of the relation between 
Rf 2 and S l (A,x,Q 2 ). These two parameterizations were used earlier to estimate the effect 
of shadowing on cc and bb production in central collisions J5| with no spatial dependence 
assumed for the shadowing. 

The Erst parameterization is a fit to recent nuclear deep-inelastic scattering data. The 
fit does not differentiate between quark, antiquark, and gluon modifications and does not in- 


b We note that at even smaller values of x. shadowing within the nucleon itself is expected III- 
However, at RHIC energies, this very low x region is not expected to be reached. 
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elude evolution in Q 2 . Therefore it is not designed to conserve baryon number or momentum. 
We define Rp 2 = Si(A,x) [[20j with 

1 + 0.0134(l/x- l/xsh) 


R 


x < x sh 


Si(A,x) = < 


l + 0.0127A 01 (l/x- l/x sh ) 

®emc b emc X Xsh X <t Xf erm i ■ (14) 

1 _ , rr . \ 0.321 
x term] \ 


R 


x 


-J 


Xfermi < X < 1 


where Rg flemc ^emc^shi Rf ®emc ^emc-^fermij ^emc 0.525(1 A ;/ 1.145v4 ! ~t~ 
0.9371 -1 + 0.88A -4 / 3 — 0.59A -5//3 ), and a emc = 1 + b emc x emc . The fit fixes x s h = 0.15, 
iemc = 0.275 and Xf erm i = 0.742. Thus, the nuclear parton densities are modified so that 


ff{x,Q 2 ) = S 1 (A,x)ff(x,Q 2 ). 


(15) 


The second parameterization, S^A, x, Q 2 ), modifies the valence and sea quark and gluon 
distributions separately and also includes Q 2 evolution [|2l|, but is based on an older fit to 


the data using the Duke-Owens parton densities |22|]. The initial scale for the evolution 
is Qo — 2 GeV and the Q 2 evolution is studied with both the standard Altarelli-Parisi 
evolution and with modifications due to gluon recombination at high density. The gluon 
recombination terms do not strongly alter the evolution. In this case, the nuclear parton 
densities are modified so that 


fy{ x i Q 2 ) 

= S?(A,x,Q 2 )F(x,Q 2 ) 

(16) 


= Sl(A,x,Q 2 )fl(x,Q 2 ) 

(17) 

f£(x, Q 2 ) 

= S$(A,x, Q 2 )/S(x, Q 2 ) , 

(18) 


where fy = uy + dy is the valence quark density and fs = 2 (u + d + s) is the total sea 
quark density. We assume that and Sf affect the up, down, and strange valence and sea 
quarks identically. The ratios were constrained by assuming that Rp 2 ~ S% at large x and 
Rf 2 ~ S2 at small x since xfy(x, Qo) — > 0 as x — ► 0. For the gluons, we take Rf 2 ~ Sy f° r 
all x JZTJ , since one might expect more shadowing for the sea quarks, generated from gluons, 
at small x. These parton densities do conserve baryon number, /d dx fy A (x,Q 2 ) = 3, and 
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momentum, /q 1 dxx(fy A (x, Q 2 ) + fs' A (x,Q 2 ) + fQ A (x,Q 2 )) = 1. at all Q 2 . We have used 
the MRS D—' and GRV HO densities with S\ instead of the Duke-Owens densities, leading 
to some small deviations in the momentum sum but the general trend is unchanged. 

Since the shadowing is likely related to the nuclear density, it should also depend on 
the spatial distribution of the partons within the nucleus so that S l (A,x,Q 2 ,r, z) —> 1 as 
s —► oo. The reduced shadowing is reasonable since the shadowing mechanism should be less 
effective when the nuclear density is low. This spatial dependence should also be normalized 
so that T f d 2 rdzp(s)S l (A, XjQ 2 ,r,z) = S l (A,x,Q 2 ) to recover the deep-inelastic scattering 
results which do not have any explicit impact parameter dependence. This approach may 
fail when x —> 1, because then the change in the structure function is likely due to Fermi 
motion, which should not exhibit spatial dependence. 

One natural parameterization of the spatial dependence follows the nuclear matter den¬ 
sity distribution, 


S' y/s = S i (A,x,Q 2 ,r,z) = l + N, 


ws 


S\A,x,Q 2 


1 + exp((s - R.A)/d ) 


(19) 


— 1 + ^wsi^iA, x, Q 2 ) — 1) 


P(s) 

Po 


where 7V W S = 1-317 is needed for the normalization to S l (A, x, Q 2 ). This form of the spatial 
dependence has a rather weak dependence on s until the nuclear surface is approached. 
Note that when s —> 0, 5'( vs < S l in the shadowing and EMC regions while Sfa s > S l in the 
antishadowing region. 

The actual spatial dependence of shadowing may be stronger if the shadowing effect is 
not directly related to the nuclear matter density distribution. This can occur if the gluons 
are not well localized within the nucleus. One can alternatively assume that the shadowing 
is related to the nuclear thickness at the collision point, proportional to the distance a parton 


from one nucleus travels through the other |23[. Therefore we also consider 


S l R (A,x,Q 2 ,r,z ) = 


1 + iV R (S‘(X X, Q 2 )- 1)^1 - (r/R A y r<R A 


( 20 ) 


r > R. 
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where Nr = 1.449 assures the normalization after the average over p(s). Similarly, when s —► 
0, < S‘ in the shadowing and EMC regions while > S' 1 in the antishadowing region. 

The normalization is higher here because of the larger region over which the suppression 
due to shadowing is reduced relative to S^s- 

We calculate the cc production cross sections in peripheral nuclear collisions with 
S l (A, x, Q 2 ) = 1, Si, and S\. As we will show, the shape of the inclusive charm quark 
Pt distributions are similar for Si and S\. Therefore, we model the spatial dependence of 
S\ only, according to eqs. (19) and (20). 


III. CORRELATION BETWEEN E r AND IMPACT PARAMETER 


Although the impact parameter is not directly measurable it can be related to direct 
observables. We discuss here the indirect measurement of the impact parameter b by means 


of the transverse energy E T |7],f24]]. Here E T = E, : y rnf + , summed over all detected 
particles in the event with masses m, and transverse momenta pxi■ It is also possible to infer 
the impact parameter by a measurement of the nuclear breakup since the beam remnants 
deposited in a zero degree calorimeter are correlated with the impact parameter. A measure 
of the total charged particle multiplicity, proportional to E T , could be used to refine the 
impact parameter determination. 

The transverse energy contains “soft” and “hard” components. The “hard” components 
arise from quark and gluon interactions above momentum po, the scale above which per¬ 
turbative QCD is assumed to be valid. Minijet production, calculated for pr ,jet > Po ~ 2 
GeV [|25| , becomes an important contribution to the dynamics of the system in high-energy 
nucleus-nucleus collisions. The hard cross section, cxff(po) = 2crj e t, twice the single LO mini¬ 
jet production cross section, can be calculated perturbatively. “Soft” processes with p T < p 0 
are not perturbatively calculable yet they produce a substantial fraction of the measured Et 
at high energies (and almost the entire Et at CERN SPS energies). These processes must 
be modeled phenomenologically. We assume a P g = u?f elastic , the inelastic pp scattering cross 
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section. Our calculation of the total Et distribution follows Ref. [^4[ . 

If the hard component is formed by independent parton-parton collisions, then the av¬ 
erage number of hard parton-parton collisions as a function of b , N AA (b), is 


= T AA (b)o%(po) . 


( 21 ) 


where cr#(p 0 ) ~ 6.5 mb at RHIC P51 and Taa(P) is the nuclear overlap function, 


T AA (b) = / d 2 rT A (r)T A (b - f) 


( 22 ) 


where the nuclear thickness function is defined as T A (r) = / dzpA(z, r). In Au+Au collisions 


at b = 0, Taa = 29/mb |26|. The Et distribution can be expressed as |2 4 


d °H _ f ,2 hy oo l N AA( b )} 

dE r J 1 N=1 N\ 


N 


— r N 

exp [-N AA (b)] / J] dE n- 

J i =i ( 


~ E ' 

a H d-t^Ti 


Ti 


(23) 


tH 


If N AA is large, dan/dE T can be approximated by the Gaussian |24 


da 


H 


dE 7 


= / d 2 b 


exp — 


AA 


(Et ~ E TH (b )) 
2a 2 H (b) 


(24) 


^Tiajjib) 

where the mean E Tl E^(b), and standard deviation, a# (5), are proportional to the first 
and second Et moments of the hard cross section, 


-fvIA 


AhW = T AA (b)a’S( Va )(E T ) 


PP 

H 


° 2 „(b) = T AA (b)a’Si^)(E 2 )’S 


(25) 

(26) 


In the rapidity interval \y\ < 0.5, a p ^(p 0 )(E t ) P h ~ 17mbGeV and a p ^ (p 0 ) (E^) p ^ ~ 


70 mb GeV 2 El 


At RHIC energies, the hard part does not dominate the soft component, proportional to 
the number of nucleon-nucleon collisions, 


N AA (b) = T AA (b)ag , 


(27) 


where a^ p ~ 30 mb. Since the soft component is almost independent of the collision energy, 


we assume, as in Ref. [gjj , that the hard and soft components are separable on the pp level 
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and thus independent of each other at fixed b. Therefore the total Et distribution is a 
convolution of the hard and soft components with total mean and standard deviation 


=AA 


Er(b) = T A A(b)[o<S(p 0 )(E T )’S + e 0 } 
o 2 (b)=T AA (b){o%(p 0 )(E$)’g + e l ] 


(28) 

(29) 


where eo and e\ are taken from lower energy data and adjusted to the same rapidity interval 
as the hard component, \y\ < 0.5, eo = 15 mb GeV, e± = 50 mb GeV 2 Shadowing, 
which affects the hard component by reducing the minijet cross section, is not included in 
these averages. Multiplying a ](f by a shadowing factor modifies the Et distribution by less 


than 10% p7|. A correction has been included here. 


Figure 2 shows the Et distribution (for y < |0.5|) for 100 GeV/nucleon Au+Au collisions 
for several different impact parameter intervals as well as the total cross section. Singling out 
a particular Et range can therefore select a rather narrow distribution of impact parameters. 
For example, requiring Et < 300 GeV selects almost exclusively events with b > R A while 
E t < 180 GeV selects events with b > 1.2 R A . 

Good event purity can be obtained with even narrower selections. For example, 300 
GeV > Et > 180 GeV largely corresponds to 1.2 R A > b > R A . An example of the 
purity can be seen in Fig. 3 which shows the range of impact parameters at Et = 200 
GeV. The distribution is centered at b = 1.27 R A with a standard deviation a ~ 0.05/)%. 
Approximately 90% of the events fall into the range 1.15R A < b < 1.35 R A , narrow enough 
to be an effective impact parameter selector. Thus at Et = 200 GeV, the impact parameter 
can be measured to within 10%. However, the statistical accuracy depends on the average 
number of collisions, proportional to Et, so that a/b ~ 1/a /Et- 

For very small Et, complications arise. The first concerns the transition from eq. (23) 

_ 

to eq. (24) which is only valid if N AA {b) is large enough for the Poisson distribution to 
be approximated by a Gaussian. For a small number of collisions, eq. (24) overestimates 
the number of low Et events, even allowing a finite probability for negative Et events. In 
practice, the agreement is quite good even at b = 1.8 R A , corresponding to T AA = 0.9/mb, 
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Naa = 5-5 and Et ~ 30 GeV. At significantly smaller Et a correction is needed. Further, 
the event by event fluctuations are large when the collision number is small, increasing the 
uncertainty in the impact parameter measurement. 

At small Et the presence of charmed quarks will alter the relationship between Et and 
impact parameter because a cc —» DD pair must have Et > 2 mo ~ 3.7 GeV. Typical values 
are Et ~ 4 — 6 GeV. Thus when Et < 20 GeV, the relationship between Et and b in charm 
events will be different. This altered relationship can be studied in simulations to correct 
the data. 

Finally, other types of interactions can contribute to charm production at low Et- The 
largest identified charm contribution in very peripheral collisions is photon-gluon fusion 


Any real detector can only measure Et in a limited rapidity interval. For example, the 
calorimeter of the STAR detector at RH1C will cover the range — 1 < y < 2 [f|J. The 
acceptance can be compensated by appropriately modifying (E t ) P h, {Et) P h, eo an d ei, given 
here for \y\ < 0.5. The accuracy scales roughly as the square root of the observed energy. A 
large acceptance can also extend the region of validity of eq. (24) to larger b. 

The non-central event selection technique to constrain the impact parameter may be 
useful in other analyses of heavy ion data. At large impact parameters, only the outer 
portions of the nuclei are involved but as the collision centrality increases, the nuclear interior 
is more deeply probed. Therefore the impact parameter variation roughly corresponds to 
the portion of the nucleus involved in the interaction, and can thus be used to study the 
difference between the parton constituents of the nuclear core and those near the surface. 


IV. RESULTS 

The best way to determine the gluon momentum fraction is to detect both charm quarks. 
Then x\ and X 2 can be fixed exactly and the shadowing mapped out. The measurements 
are relatively easy to interpret if y = —y since X\ = X 2 - After first discussing the general 
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results when the kinematic variables are integrated over, we show the px distributions for 
the MRS D—' and GRV HO parton densities assuming both the c and c are detected. The 
low experimental efficiency for detecting charm suggests that it is unlikely for both quarks 
to be detected in an event. Thus we subsequently discuss the feasibility of the study if only 
one of the charm quarks is detected. 

Figure 4 shows the cc production cross section as a function of impact parameter for 
b > Ra with S = 1, S\ and S 2 at RHIC f|. The cross sections were calculated by integrating 
eq. (1) over the c and c four-momenta. The rates for these non-central collisions are still 
quite large. Without shadowing, for b > 1.2 Ra the charm cross section is 2.9 b while for 
b > 1.8 Ra it is still 200 mb. At the RHIC Au+Au design luminosity, 2 x 10 26 cm -2 sec -1 0, 
this results in 6300 and 430 million cc pairs/year (3000 hours). Thus these measurements 
will not be statistics limited, even with the roughly 35% reduction in cross section when 
shadowing is included. 

Figures 5 and 6 show the charm quark px distributions in two different impact parameter 
intervals, b > 1.2 Ra, roughly corresponding to E T < 180 GeV in Fig. 2, and b > 1.8 Ra, 
for several selected c and c quark rapidities. The results with the MRS D—' and GRV 
HO parton densities are compared. By measuring charm quarks as a function of px for 
a variety of rapidities, different values of x\ and X 2 are probed. For example, px = 0, 
y = y = 0 corresponds to x\ = x 2 = 1.3 x 1CT 2 while p T = 0, y = 2 and y = — 2 corresponds 
to X\ = x 2 = 5.1 x ICR 2 . At px ~ 2.1 GeV X\ and x 2 are doubled, moving into the 
antishadowing region for \y\ = 2. Thus varying X\ and X 2 changes the relative strength of 
the shadowing. Calculations with S = 1, Sf, S qws, Si,r and S 2 are shown in each case. 

In every case considered, the unshadowed cross section is larger than the shadowed cross 
sections. The total cc production cross sections with Q oc nm c differ only by 2% in pp 
collisions. (Recall that n — 2 for MRSD—' and n — 1 for GRV HO.) When the total cross 
section is computed by integrating an inclusive cross section where Q oc nmx, the difference 
increases to ~ 6% due to the running scale in the parton distributions and a s . The inclusive 
distributions reflect the low x and Q 2 behavior of the parton distributions. The MRS D—' 
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gluon distributions are always decreasing as a function of px- However, the GRV HO gluon 
distributions are still valence-like at low Q. Thus for y = y = 0 and px < 1-5 GeV the gluon 
distribution continues to increase, causing the observed ~ 15% difference between the S = 1 
distributions at pt ~ 0 in Figs. 5(a) and (d). At larger rapidity and x, such as in Figs. 5(c) 
and (f), the difference is reduced to ~ 8%. 

The shadowing functions affect the charm px distributions differently for the MRS D—' 
and GRV HO parton distributions because of the difference in the scale Q 2 . In general 
increases more rapidly with x than S\ between the shadowing and antishadowing regions. 
With the MRS D—' parton distributions, at px ~ 0, Si ~ S'f for Q fa 2 m c . As px increases, 
S 2 > S 1 due to the evolution of .SV Therefore when px ~ 1 GeV, the px distribution with 
S 2 will be ~ 10% larger than the distribution with S\. This continues to hold as px rises, 
as shown in Figs. 5(a), (b) and (c). The GRV HO case is different because of the lower 
scale. There, the evolution of S 2 with Q 2 does not begin until Qo = 2 GeV, corresponding 
to px ~ 1.5 GeV. For p T < 1.5 GeV, Si > Stf. At p T ~ 2 GeV for y = y = 0, the evolution 
of S 2 causes the situation to be reversed and > Si, as can be seen by inspection of Fig. 
5(d). At larger rapidities, the larger slope of S'f in the shadowing region cause the switch 
between S± and S 2 dominance to occur at lower values of px, even before the evolution of 
S 2 begins, since x is larger at small px and large y. 

Including spatial dependence in S\ increases the cross section toward the S = 1 value at 
high b where the nuclear density is low. The cross section is now larger because the lower 
density near the nuclear surface reduces the shadowing. As the impact parameter rises, 
the tails of the density distributions are probed and the shadowed cross sections approach 
the S = 1 result. This happens relatively slowly, especially for S gws, since the density 
is nearly constant except within d of the surface. The shadowing is thus almost constant 
except near the nuclear surface. For gold, d = 0.535 fm while b = 1.8 Ra, the lower bound 
on the impact parameter in Fig. 6, corresponds to collisions within 1.2 fm of the surface so 
that some collisions occur below the surface layer in at least one nucleus. In both Figs. 5 
and 6, > Agws because the dependence on the nuclear thickness (albeit for a spherical 
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nucleus) decreases the effects of shadowing already at small r while Syws is almost constant. 
The effect is more apparent for larger impact parameters. When b > 1.2 Ra, both spatial 
forms increase the cross section about 15% over S\. For b > 2 Ra the spatial results are 
approximately halfway between the cross sections with 5 = 1 and 5 — S\. The similarity 
of results between the two spatial parameterizations suggests that the parton localization 
measurements may not be too hard to interpret. 

Thus measurements of charm quark production at large impact parameters probe the 
nuclear surface where shadowing effects are greatly reduced, and, for extremely peripheral 
collisions, the limit of independent pp collisions is regained. As the collisions become more 
central, the charm quark production rate should begin to deviate from the naive expectation 
from superimposed pp collisions. By measuring charm production as a function of impact 
parameter, it is possible to watch the shadowing turn on with the rate of increase providing 
a measure of parton localization in the nucleus. 

So far we have assumed that both the c and c quarks are detected. Given the low 
efficiency for detecting charm quarks, either by their semileptonic decays or by reconstruction 
of specific final states, it is worth considering what can be learned if only one of the quarks 
is detected. Fig. 7 shows the rapidity distribution of the c quark, assuming that the c quark 
is detected at y — 0 and pt = 0 assuming 5 = 1, S±, S 2 , 5i t ws and 5 i ; r. Kinematically, this 
situation corresponds to charm pair invariant mass M 2 = 2m 2 (1 +coshy) so that increasing 
y corresponds to increasing phase space along with increasing invariant mass. The cross 
section increases until y ~ ±1 where M « 3.4 GeV and decreases with larger M, typical 
for invariant mass distributions JT2|. Fig. 8 shows the single charm p t distribution at y = 0 
integrated over y for b > 1.8 Ra- The results are similar to the case when both quarks are 
detected. Although some information is lost if only a single quark is detected, the trends 
remain the same as those seen in Fig. 6. Therefore it should still be possible to extract the 
shadowing information from the data. 
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V. DISCUSSION 


If the charmed quark rapidity and momentum can be measured over a broad range of 
impact parameters, the gluon momentum distribution and its spatial/density dependence 
can be measured. However, there are a number of difficulties involved in relating these 
calculations to measurements. Charm is normally detected either via its scmileptonic decays 
or through reconstruction of selected decay modes. While the detection of leptons from 
semileptonic decays is fairly straightforward, the lepton px and y differ from that of the 
parent hadron. The parent hadron distribution can also differ slightly from that of the 
initially produced quark although the hadronic environment reduces this effect J3(|. While 
this momentum shift does not create any fundamental problems, it adds another intermediate 
step which must be correctly modelled. Fully reconstructed charm decays such as D* + —> 
D°tt + —> (K~ 7r + ) 7r + could allow for a full reconstruction of the meson direction, reducing 
the uncertainty in the determination of the charmed quark px and y. However, the small 
branching ratios and low efficiency for detecting these decays probably preclude the useful 
detection of both charmed quarks in a pair. 

In addition to gold, RHIC will accelerate a variety of lighter nuclei. The surface layer is 
a larger fraction of the nuclear radius in lighter nuclei. In this case, the Woods-Saxon and 
square root spatial dependencies should more closely match over the full range of impact 
parameters. Since RHIC is also a pA collider, the gluon localization could in principle 
be probed for an individual nucleus. However, for pA, the number of collisions is small 
enough for the Gaussian approximation to break down, rendering the Ex to b correlation 
problematic. The A dependence of charm production at various impact parameters can 
in any case provide an additional handle on interplay between shadowing and its spatial 
dependence. For pA, dileptons can also be used to probe gluon shadowing 0. 


At LHC, similar calculations can be made for cc and bb production. The higher energy 
implies that the charm and bottom pairs will be produced at much lower x, increasing 
the importance of shadowing and further reducing the production cross sections. Thus the 
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sensitivity of the cross section to the spatial dependence will be enhanced. 

VI. CONCLUSIONS 

We have calculated charmed quark production in non-central Au+Au collisions for several 
different structure functions and assumptions about nuclear shadowing. 

Shadowing reduces the charm production cross section up to 35%. However, when the 
spatial dependence of shadowing is included, the effect is decreased. By measuring the 
charmed quark production rates as a function of impact parameter, it is possible to study 
the effect of shadowing and its localization within the nucleus. This spatial dependence 
provides an indication of the gluon recombination distance scale. 

The correlation between impact parameter and transverse energy has been used to fix 
b. We have shown that the impact parameter determination is reliable to within a 10% 
statistical uncertainty on an event-by-event basis for b « 1.2 R&. 
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FIG. 1. The collision geometry of nuclear collisions in the plane transverse to the beam. The 
parton-parton collision point is indicated by A and b is the impact parameter. 
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FIG. 3. Distribution of impact parameter for events with Et = 200 GeV. 
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FIG. 5. The pt distribution of cc pairs for the MRS D—' (a), (b) and (c) and GRV HO (d), (e) 
and (f) parton densities. We select events with b > 1.2 Ra and 3 sets of c and c quark rapidities: 
y = 0, y = 0 in (a) and (d); y = 0, y = 2 in (b) and (e); y = 2, y = —2 in (c) and (f). The solid 
curves are for S' = 1. The spatially independent shadowing results are given by the dashed (Si) 
and dotted (S 2 ) curves. The effect of the spatial dependence on Si is also shown. The dashed 
curve with the filled squares shows the result with SgR and the dashed curve with the open circles 
gives the result with Sgws- In (a), (b), (d) and (e) the Si and Srivs curves overlap. 
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FIG. 7. The c rapidity distribution for pr = 0 and the charm quark is produced at y = 0. The 
solid curve is with 5 = 1. The spatially independent shadowing results are given in the dashed and 
dotted curves for S\ and S 2 respectively. The effect of the spatial dependence on 5 1 is also shown. 
The filled squares shows the result with 5 i ; r and the open circles gives the result with 5i ; ws- 


27 







FIG. 8 . The pr distribution for single charm quarks with y = 0 for MRS D—' (a) and GRV HO 
(b) parton densities. We have selected events with b > I. 8 -R 4 . The solid curves are with 5 = 1. 
The spatially independent shadowing results are given in the dashed and dotted curves for S\ and 
S 2 respectively. The effect of the spatial dependence on S\ is also shown. The dashed curve with 
the filled squares shows the result with Srr and the dashed curve with the open circles gives the 
result with S^ws- 





